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We discuss transport properties of dynamical fluid composed of quasi-particles whose masses 
depend on temperature and charge chemical potentials. Based on the relativistic kinetic theory 
formulated under the relaxation time approximation, we derive a general expression for the bulk 
viscosity in the quasi-particle medium. We show that dynamically generated particle masses imply 
an essential modification of the fluid compressibility. As an application of our results we consider a 
class of quasi-particle models with the chiral phase transition belonging to 0(4) and Z(2) universality 
class. Based on the Ginzburg-Landau and the scaling theory we study the critical properties of the 
bulk viscosity £ near the phase transition. We show that extrapolating the results of kinetic theory 
under the relaxation time approximation near the critical region the £ does not show singular 
behavior near the 0(4) and Z{2) critical point through static critical exponents. 

PACS numbers: 25.75.Nq, 24.85.+p, 12.39.Fe 
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1. INTRODUCTION 

In the hydrodynamical evolution of fluid leading dissi- 
pative processes can be quantified by the transport coef- 
ficients, the shear r\ and bulk £ viscosities. Their values 
and properties are not only carrying information on how 
far the system appears from an ideal hydrodynamics but 
can also provide relevant insight into the fluid dynamics 
and its critical phenomena [H, H, S IE S IzL HI • For cer- 
tain materials, e.g. helium, nitrogen or water, the shear 
viscosity to entropy ratio rj/s is known experimentally to 
show a minimum at the phase transition |3[ . On the other 
hand the bulk viscosity (/s was argued to be large or even 
divergent at the critical point [E [U, @, Q ■ The recent 
Lattice Gauge Theory calculations seem to be consistent 
with the expectation of decreasing rj/s and increasing 
C/s toward the QCD phase transition from above fflfiol]. 
Thus, transport coefficients are of particular interest to 
quantify the properties of strongly interacting relativistic 
fluid and its phase transition [3(. 

In modeling strongly interacting media near equilib- 
rium, the interactions usually lead to a quasi-particle de- 
scription with its mass depending on thermal parameters 
fill ]. The thermodynamics of perturbative QCD is also 
to large extent quantified by thermal quark and gluon 
masses which are T and ji dependent [12||. In the QCD- 
like chiral models where the phase transition is governed 
by an effective mass generated through the dynamics, the 
quasi-particle mass plays a role of an order parameter and 
thus is sensitive to change in thermal parameters. 

In this paper we discuss the bulk and shear viscosities 
of the dynamical fluid composed of quasi-particles with 
T and fi dependent masses. Our calculations arc based 
on the kinetic theory in the relaxation time approxima- 
tion. We derive a general expression for the viscosities in 
the quasi-particle medium which has a broad spectrum 
of applications. At vanishing chemical potential or for 
fixed particle masses our results are consistent with that 



formulated in |l3| and [14( respectively. 

Extrapolating the calculated bulk viscosity to the 
phase transition region we discuss its critical properties. 
Based on the Ginzburg-Landau and the scaling theory we 
show that under the relaxation time approximation the 
£ is not expected to show a singular behavior near the 
0(4) and Z(2) critical point through the static critical 
exponents. 



2. SHEAR AND BULK VISCOSITIES FROM 
TRANSPORT THEORY 

The transport parameters, the shear 77 and the bulk £ 
viscosities, are defined as coefficients of the space-space 
component of a deviation of the energy momentum ten- 
sor from equilibrium. In a medium composed of bosons 
and/or fermions with the momentum distribution func- 
tion f(p, x) for a particle and /(p, x) for an anti-particle 
the energy momentum tensor is defined as 



(IT 



E 



[f + fl 



(2.1) 



where dT = gd 3 p/ (2ir) 3 is the integration measure in the 
momentum space with the degeneracy factor g associated 
with the particle quantum numbers. The four momen- 
tum = (E,p) with E(p) = yjp 1 + M 2 and M being a 
particle mass.^ 1 . 

Assuming that the medium appears near equilibrium 
we introduce the particle momentum distribution /(p, x) 
in the following form 



f = (e 



(E-p-u=FM)/T 



±1) 



(2.2) 



# x In general, it is not necessary to specify this dispersion relation. 
The following derivation is valid for any E(p) 



2 



where u is the flow velocity and p is the chemical po- 
tential related with any conserved charges. The ±1 cor- 
responds to fermion and boson statistics whereas =F/i to 
particle and antiparticle contributions, respectively. 

In our further discussion we assume that M is not 
necessarily a bare particle mass but rather a dynamical 
quasi-particle mass which depends on temperature and 
chemical potentials, thus M = M(T,p) in Eq. (|2.2[) . 

A deviation of the system from equilibrium is quan- 
tified by the corresponding change of distribution func- 
tions Sf = f — /o with /o being the equilibrium particle 
momentum distribution. In the relaxation time approxi- 
mation [l5j , the Sf is obtained from 



P^fo = -—Sf, 



(2.3) 



where the collision time 



The pressure P is a function of T and p, thus 

dP _ dP dT dP dp 
7ft ~ df~dt + ~dp~ 7ft ' 



(2.11) 



therefore one arrives at the following equations for the 
time dependence of T and p; 



or T fdP\ 



dp 

7ft 



dP 
7ft 



dP 
dn 



V-u. (2.12) 



Consequently, the change of the energy momentum ten- 
sor (|2.5[) becomes 



nf(T, p)(va(T, p)) 



(2.4) 



is determined by the thermal-averaged total scattering 
cross section (va) , with the relative velocity of two collid- 
ing particles v and the particle density nf in equilibrium. 
The Sf results in the corresponding change in the energy 
momentum tensor 



5T^ = 6T? V + ST$ V , 



f 



drr^f (l±f )q f (p;T,fx), (2.13) 



where we have introduced 



ST^ V 



dr^rP a d a 



rf + f/o 



(2.5) 



The thermodynamic quantities arc time dependent 
through thermal parameters. The time dependence of 
T and p are obtained from the charge number density 
doj° = and the energy density doT 00 = conserva- 
tions. The charged current = is defined by 



[/-/] , 



The energy and charge density conservations can be ex- 
pressed in terms of the thermodynamic quantities as 



de 
7ft 



dn 
7ft 



(e + PW -u 



= -IT 



dP dP 



dT 
nV • u = 



V • u. 



dp 



(2.7) 
(2.8) 



where the energy e and charge n densities are related 
with the pressure P through the thermodynamic rela- 
tions: e = TdP I dT — P + pdP/dp. From the chain rules 



dP 
~dt 



dP de 
7ft 7ft 



dP dn 
dn dt 



(2.9) 



P 1 , fdP\ ( T dE dE 



' . dn ) £ ( dp, T 1 



dkUl S kl -^W kl , (2.14) 



with the following tensor decomposition 



(2.6) 9 k u l = i (d k u l + dm k - ^S M diU^j + ^SudiU* 



\w kl + y k Au i . 



(2.15) 



With the above decomposition, the change of the tensor 
T 1 - 7 can be written as a sum of traceless Wa and scalar 



part 



61*' = -(S tJ d k u k - r)Wij , 



(2.16) 



which defines the transport coefficients, the bulk £ and 
the shear n viscosities respectively. Implementing in 
Eq. (|2.13[) the energy conservation 

J dTE [r/ (l ± f )q f + r/ (l ± h)qj\ = (2.17) 



and from Eq 
~dP 



one gets 



dP 
7ft 



dP 



T- 



de \ dT 



dP 
dp 



dP dP 
dn dp 



V-u. (2.10) 



leads to the final expression for the transport coefficients 
obtained under the relaxation time approximation in the 
quasi-particle model with T and p dependent masses. 
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The bulk £ = C/ + C/ an d the shear r\ = T}f + r}j viscosities in a medium composed of one type of particle/ antiparticle 
is obtained as 



1 f d 3 p p 4 
15T J (2tt) 3 E 2 



[ 5 r/o(l±/o)+5r/ (l±/o)] , 



(2.18) 



and 



d 3 p 



A/ 2 



(sr/o(l±/o)+St/o(l±/o)) 

M 2 - - f ftP 

— ( 5 r/ (l ± f ) - ,9f/ (l ± /„)) ( — 



/jr^ _ (ftP_ 
\3E ~ \ fte 



rp T 0E 8E 

E - T ftf-^ 



ftP\ ftE 
ftn ) ft [i 



(2.19) 



respectively. The above results can be obviously generalized to any systems composed of different particle species by 
summing up their contributions in Eqs. (|2.18jl and (|2.19j) . 



The derivatives of pressure dP/de\ n and ftP/dn\ e en- 
tering in Eq. (|2 . 1 9[) can be expressed in terms of the 
net particle number density n, the entropy density s and 
different susceptibilities Xxy which are defined by 



ftP 
d/j, 



ftP 
df ' 



Xxy 



d 2 p 

ftxfty 



(2.20) 



Applying the Jacobian methods to the above derivatives 
of P one finds 



ftP 

ftP 
ftn 



C'vXw 

nTxTT + {nn - sT)x^t - snxw 



(2.21) 



CvXw 



(2.22) 



with Cy being the specific heat at constant volume and at 
constant s/n. The Cy can be expressed through different 
susceptibilities as 



T 



ds_ 
df 



The shear viscosit 
viously obtained in 



n 



T 



Xtt 



xIt 



(2.23) 



1G] 



has the same form as pre- 
where the thermal modifi- 



cation of particle dispersion relations was not included. 
Thus, the shear viscosity is not directly affected by the 
quasi-particle dynamics. However, the bulk viscosity is 
essentially modified by the terms ftE/ftT and dE/d/j, that 
appear only when there is an explicit (T, /^-dependence 
of the particle mass. For M(T,fj,) = const., Eq. (|2.19p 
is reduced to the result obtained in 

0. At vanishing 

chemical potential Eq. (12.19J coincides with the expres- 
sion recently formulated in [l3j for an interacting quark- 
gluon plasma. 



3. TRANSPORT COEFFICIENTS NEAR THE 
PHASE TRANSITION 

To quantify the transport properties of thermodynamic 
systems characterized by viscosity coefficients one would 



need to formulate a specific model for particle interac- 
tions. However, there arc some generic properties of 
£ which can be discussed in a model-independent way 
through the universality arguments. In this context, of 
particular interest is the specific behavior of £ near the 
phase transition. For the bare particle mass, the dy- 
namics is entering only through the collision time. Thus, 
information on phase transition can only be contained in 
the change of the thermal averaged cross section. How- 
ever, in the quasi-particle picture due to dynamical par- 
ticle masses, further information on the phase transition 
appears in £ through an explicit contribution of differ- 
ent observables which are sensitive to critical phenom- 
ena. From Eq. (|2.19[) it is clear that such observables 
are the specific heat CV and different susceptibilities. In 
addition, in the model where the particle mass M(T, /z) 
is related with an order parameter, the thermal deriva- 
tives of M can be as well singular at the phase transition. 
This is e.g. the case when considering effective models 
related with the chiral phase transition 0, [TtJ ■ 

To study the sensitivity of the bulk viscosity to the 
phase transition we separate from Eq. (|2.19p the term 
C (dcr) that can be singular near the critical point, 



/•(dcr) _ _ J_ f 

3T J (2tt 



d 3 p M 2 



ftP 



f E 
dE 



E-T 



[( ff r/o(l±/o)+5f/ (l±/o)) 
dE\ fftP\ ftE" 



ftT 



/' 



9/x 



dP 
ftn 



/ftP 

-( 5 T/o(l±/o)-fff/ (l±/o)) 



(3.1) 



Hereafter, assuming a validity of the relaxation time 
approximation we discuss the critical behavior of C' der ) 
in different models with a phase transition. In particular, 
based on the Ginzburg-Landau as well as on the scaling 
theory we analyze the scaling of £ in the 0(4) and Z(2) 
universality classes which are expected in the QCD chiral 
phase transition. 
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3.1. Mean field scaling of the bulk viscosity in the 
Ginzburg-Landau model 

According to the Ginzburg-Landau theory, close to the 
phase boundary the thermodynamic potential may be ex- 
panded in a power series of the order parameter M, which 
plays the role of the dynamical particle mass 

fi(T > /i)~n (T, W M = 0) 

+ M 2 + K£ltl M 4 + £ M 6 _ mM (3>2) 

with c > 0. The order parameter is determined in such 
a way that the free energy is minimized. Thus, M is 
obtained as the solution of the gap equation for m = 



M 2 = — ± — Vb 2 -Aac. 
2c 2c 



(3.3) 



The fluctuations of M arc defined through the suscepti- 
bility 



X 



dM 
dm 



3bM 2 + 5cM 4 



(3.4) 



With the particular choice of the parameters, a = for 
b > the thermodynamic potential describes the second 
order phase transition. On the other hand for a = b = 
the system exhibits the tricritical point (TCP). Thus, the 
Ginzburg-Landau model has generic critical properties 
expected in the QCD chiral phase transition. 

In the vicinity of the phase transition the tempera- 
ture and chemical potential dependence of the coefficients 
a(T, jj) and b(T, Li) can be parameterized as 



a(T, fx) = a\T - T c \ + /3\fx - fx c 



(3.5) 



with constant a and j3. 

From the thermodynamic potential (|3.2[) one gets all 
relevant thermodynamic quantities near the phase tran- 
sition. In particular, near the TCP the singular part of 
susceptibilities 



Xtt 



op 

T ' XaiT ~ ~b ' 

and derivatives of the dynamical mass 

dM I a dM 1/3 
~d~T M~b' ~d/7 Mb 



(3.6) 



(3.7) 



With the above scaling relations and from Eqs. (|2.21|) 
and (|2.23p one finds, that near the TCP the pressure 
derivatives and the specific heat are finite and that the 
potentially singular part £( dcr ) behaves as 



c 



Consequently, due to the vanishing coefficient of the sec- 
ond bracket in Eq. (|3.8p together with the scaling [13] 



AP 
~b~ 



M\ ~ t c 



(3.9) 



the singularities of the susceptibilities do not show up in 
C near TCP. 

The singular part of £( dcr ) also vanishes at the second 
order transition since there C (der) ~ M 2 — > 0. In addi- 
tion, if the explicit chiral symmetry breaking is absent 
then also the regular part of ( vanishes. Thus the £ is 
precisely zero at the critical point. 

The above example shows that within the mean field 
dynamics the bulk viscosity is non-singular at the second 
order phase transition and at the TCP. 



3.2. Scaling of the bulk viscosity in the 0(4) and 
Z(2) universality class 

In order to verify the singular behavior of £ along the 
second order line and at the critical end point (CEP) that 
belongs to the 0(4) and Z(2) universality class respec- 
tively, one needs to go beyond the mean field approxi- 
mation. The scaling behavior of different observables in 
the vicinity of phase transition emerges from the scaling 
function of the singular part of the free energy. 

Along the 0(4) transition line and for \iJT << 1, the 
singular part of the free energy can be parameterized as 

m 



where 



F s (T,fi)~t 2 - a f s (l,t-P 5 h) 



t = t + Ap 2 



(3.10) 



(3.11) 



with t = \T — T c \/T c , ft = n/T c and h being an external 
field. The a, (3 and S are the critical exponents in the 
0(4) universality class. 

From the scaling function (|3 . 1 Of) . by taking derivatives, 
one finds the properties of all relevant quantities required 
in Eq. (|2. 19|) to find the behavior of the bulk viscosity 
near the 0(4) line. In the following we consider only the 
transition point at fi = 0. From Eq. (|3 . 10[) one gets 



Xtt ~ t 



r 



M ~ , 8M/8T - t^- 1 



(3.12) 



Substituting the above scaling to Eq. (|2.19|) one finds 
that near the 0(4) transition point the singular part of 
the bulk viscosity £( dor ) scales as 



a+4/3-1 



(3.13) 



Thus, with the 0(4) critical exponents: a ~ —0.24 and 
(3 ~ 0.38, the singular part of ( vanishes at the critical 
point. Consequently, in this approach, there is no singu- 
larity in the bulk viscosity along the 0(4) line. 

The above scaling of C can be different near the CEP, 
because it belongs to the Z(2) universality class of the 
three dimensional Ising model. To match the spin system 
parameters, the reduced temperature t and the external 
field h, to those in the QCD chiral models we follow the 
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discussion of [19( and replace: t — > a t t + btp, and h — > 
a/jf + bhp,. The singular part of the free energy in the 
Z(2) universality class is parameterized as: 



F s {t,»)~hr¥f s {h- l /e 5 t,i). 



(3.14) 



This scaling form of the free energy leads to the following 
critical behavior of the thermodynamic quantities near 
CEP at t = 0: 



XjU/i.TT./iT 



h -i/ps _ 
dM/dT 



C v ~ /i 



-7//3«5 



(3.15) 



Consequently, the singular part of the bulk viscosity 
()2.19[) in the Z(2) universality class scales as 



(3.16) 



Substituting into Eq. (|3. 16|) the Z(2) exponents: (3 ~ 
0.31, <5 ~ 5.2 and 7 ~ 1.25 one finds C ~ h 54 . Thus, 
similarly as in the 0(4) universality class the bulk vis- 
cosity is non-singular near the CEP. 

The above scalings of £, obtained under relaxation time 
approximation, are different from that recently found 
in Q. There, it was argued that the critical behavior 
of the bulk viscosity is governed by the critical exponent 
of the specific heat since £ ~ Cy # 2 . Consequently, the 
bulk viscosity was shown to diverge at the CEP in the 
Z(2) universality class. In our approach, under the relax- 
ation time approximation, the scaling of £ is determined 
by the product 



(der) M^fdM dM\ 
^ ~ C v \dT dfi J 



(3.17) 



Consequently, C is proportional to Cy 1 rather than to 
Cy . The same, negative power of Cy is also reported 
in [2(|. The bulk viscosity contains both the regular 
and derivative terms: the former is alway positive and 
the latter can change its sign dependently on the model. 
However, the overall sign of ( should be always posi- 
tive for any T and ji. This condition could be used as 
a constraint and a limitation on the thermal parameter 
dependence of dynamical quasi-particle masses to satisfy 
the H-theorem. 

From the above discussion, valid under relaxation time 
approximation, one sees that £ is not expected to show 
divergent behavior near the critical end point through 
the static critical exponents. However, an extension of 
the relaxation time approximation, that allows to incor- 
porate properly a long-range fluctuations of soft modes, 
can result in divergence of ( at the CEP as well as at 



the second order 0(4) transition (20[. In this case, diver- 
gence of the bulk and shear viscosities is governed by the 
dynamic rather than static critical exponents [111, [22] ^ 3 . 

In the presence of the dynamical critical phenomena 
the bulk viscosity was shown to scale as ( ~ r z " +a at the 
second order phase transition, with the dynamical critical 
exponent z, determining the critical slowing down of the 
system relaxation time t ~ t~ zv and with v and a being 
the static critical exponents of the correlation length and 
the specific heat respectively [22]. Consequently, when 
approaching towards the Z(2) critical point the bulk vis- 
cosity diverges due to the large value of z ~ 3 expected 
in the dynamical universality class of the 3-dimcnsional 
Ising model [21,23]. The ratio, £/t ~ Cy 1 , however, is 
proportional to the inverse of the specific heat and thus 
vanishes at the critical point. This scaling property of the 
ratio is consistent with our finding in Eq. (40) obtained 
within the relaxation time approximation. However, due 
to the fact that in the relaxation time approximation 
the r stays finite even at the 2nd-order critical point the 
bulk viscosity obtained in our approach does not include 
a dynamical critical behavior which is quantified by the 
dynamical critical exponent. 



4. CONCLUSIONS 

We have studied the non-equilibrium properties of a 
quasi-particle medium at finite temperature and density 
using the kinetic theory in the relaxation time approx- 
imation. Assuming, that the quasi-particle masses are 
temperature and chemical potential dependent, we have 
derived a consistent expression for the bulk viscosity co- 
efficient f. We have shown that in the presence of dy- 
namical mass M the fluid compressibility is essentially 
modified. Our result for £ is valid in different physi- 
cal systems where interactions yield some modification 
in particle dispersion relations through explicit variation 
of M with thermal parameters. 

We have applied our results to a class of effective chi- 
ral models where the dynamical mass is identified as the 
order parameter for the chiral phase transition. We ex- 
trapolated the bulk viscosity obtained in the relaxation 
time approximation to the phase transition and studied 
the influence of critical fluctuations on £. We have shown 
that under the mean field dynamics as well as in the pres- 
ence of quantum fluctuations implemented through the 
scaling functions, the bulk viscosity is not sensitive to the 
chiral phase transition. This is because, the singularities 
generated from the susceptibilities arc totally canceled 
in £ at the 0(4) as well as at the Z{2) critical point. A 
possible modification of our conclusions due to a proper 
treatment of soft modes at the phase transition is not ex- 



# 2 For the recent discussion of a possible problem with such scaling * 3 The dynamic universality class of the QCD critical point is that 
behavior of the bulk viscosity see e.g. Ref. of the model H [23ll . 
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eluded. Nevertheless, the scaling properties given in this 
paper indicate a tendency of the transport coefficients 
when approaching the critical point. The phenomeno- 
logical relevance of soft modes at the phase transition 
strongly depends on how large the critical region is. If the 
critical region where the description through the static 
critical exponents is not valid anymore is very narrow in 
thermal parameters, it may be rather hard to observe the 
singularity of bulk viscosity around the critical point. 
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